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Fully resolved simulation data of flow separation over 2-D humps has been used to 
analyze the modeling terms in second-moment closures of the Reynolds-averaged Navier- 
Stokes equations. Existing models for the pressure-strain and dissipation terms have been 
analyzed using a priori calculations. All pressure-strain models are incorrect in the high- 
strain region near separation, although a better match is observed downstream, well into the 
separated-flow region. Near-wall inhomogeneity causes pressure-strain models to predict 
incorrect signs for the normal components close to the wall. In a posteriori computations, 
full Reynolds stress and explicit algebraic Reynolds stress models predict the separation 
point with varying degrees of success. However, as with one- and two-equation models, the 
separation bubble size is invariably over-predicted. 
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Reynolds stress anisotropy, 2uiUj/uiUi — 25ij/3 
flatness of stress anisotropy, 1 — 9/8[II a — III a ] 
pressure diffusion 
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turbulent kinetic energy, m 1 2 /s 2 
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turbulence production, — \upukiduj / dxk) + ujuk{dui/dxk)\ 
ratio of traces of production and dissipation tensor 
Reynolds number, Ui n h/v 
turbulence Reynolds number, k 2 /{ue) 

nondimensional mean flow strain rate tensor, k / e{dUi / dx j + dUj/d x/)/2 
turbulent diffusion 

non-dimensional, symmetric, deviatoric tensor 

maximum inlet velocity, m/s 

velocity vector, m/s 

specific Reynolds stresses, m 2 /s 2 

friction velocity, m/s 

nondimensional mean flow vorticity tensor, k/e{dUi/dXj — dU j/dx/)/ 2 
separation coordinate, 100(x — Xs)/{Xr — Xs) 

^-coordinate of start of separation 
^-coordinate of reattachment 
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dissipation (e.jj/2), m 2 /s 3 
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dissipation tensor 

pressure-strain 

kinematic viscosity, m 2 /s 

second-invariant of stress anisotropy tensor, 

second- invariant of dissipation tensor, dikdki 

third-invariant of stress anisotropy tensor, ajkCLkidij 

third-invariant of dissipation tensor, djkdkidij 


Subscript 

i,j,k,m,l,p,q indices of tensors 


Superscript 

+ wall units, scaled by inlet u T and v 

(•) Reynolds-averaged quantity 


I. Introduction 

Several years ago, widely used Reynolds-averaged Navier-Stokes (RANS) turbulence models were bench- 
marked for a 2-D hump separated flow in a workshop. 1 All RANS models — including one- and two-equation 
linear eddy-viscosity models — were inadequate in predicting the separation bubble size. In particular, when 
separation location was predicted well, reattachment was predicted too far downstream. Non-linear eddy 
viscosity models — including explicit algebraic Reynolds stress models (EARSM) — were also poor. Second- 
moment Reynolds stress models (RSM) were not tested in the workshop, but have been shown in other sep- 
arated flow studies to be similarly deficient. 2 Failure of RSMs- both algebraic and differential transport — is 
believed to be primarily caused by the models used for the pressure-strain and dissipation tensors. With 
the recent availability of highly-resolved simulation data sets (Table 1) for separation on 2-D smooth-body 
configurations, the efficacy of existing models can be studied using a priori and full computational analysis. 

Turbulent flow separation from smooth surfaces is characterized by deceleration of flow prior to flow 
reversal, followed by a separated shear layer away from the wall. The shear layer creates turbulence pro- 
duction that drives the flow away from equilibrium. RSMs provide transport equations for each component 
of the Reynolds stresses, hence providing, in theory, a higher fidelity prediction of mean flow and Reynolds 
stresses compared to 2-equation eddy- viscosity models. However, in practice they are less robust and occa- 
sionally do not agree with data as well as simpler eddy-viscosity models. The roots of the latter deficiency 
lie in the modeled terms of the Reynolds stress transport equation. The purpose of this study is to analyze 
existing RSMs in the context of 2-D separated flows, as a first step toward identifying suitable models for 
pressure-strain and dissipation terms. 

To elucidate the contributions of the RSM modeling components, consider the fully-developed channel 
flow in Fig. 1. The production ( P t j ) and viscous diffusion terms can be evaluated exactly. Four 

terms need modeling (highlighted in color in the figure); however, most modelers represent them using three 
models. Model 1 (red) accounts for turbulent-diffusion and pressure-diffusion (X) ; + D p i j), model 2 (green) 
for pressure-strain correlation ($^), and model 3 (blue) for dissipation (e^). In the log-layer region of the 
flow, the balance is primarily between the pressure-strain and production. Hence pressure-strain modeling is 
considered an important component of RSM. However, in the viscous sublayer, dissipation balances viscous, 
pressure diffusion, and pressure-strain terms. Hence all three modeling components are important, implying 
that in order to predict the stress accurately all three models need to be accurate in the near- wall region. The 
near-wall region is also inhomogeneous and poses additional difficulties to modelers; inhomogeneity effects 
are often incorporated into homogeneous pressure-diffusion and dissipation terms. A review of different 
approaches is given in Gerolymos et al. 3 

Fig. 2 shows DNS results for separated flow over a hill in a convergent-divergent channel. 4 This flow 
demonstrates some of the characteristics of 2-D separated flows. In particular, the flow is in turbulence non- 
equilibrium because the V /s in the separated shear layer is significantly greater than one. The bubble also 
encloses two regions of negative production close to the wall. Investigations by Bentaleb et al. 5 show that 
the ratio of turbulence length scale to mean-flow length scale reaches levels as high as ten in the recirculation 
zone near the wall. This ratio is a measure of the non-locality of the flow, implying that single-point closure 
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models may not be well-suited to predicting the separated flow region. Although the current work deals only 
with single-point closure modeling, this result from Bentaleb et al. suggests that further work is needed to 
investigate the inclusion non-local effects in RANS models. 

The objective of this study is to analyze the turbulence dynamics of smooth-body separated flows using 
existing Large-Eddy Simulation (LES) and Direct Numerical Simulation (DNS) data sets in an attempt to 
identify weaknesses of existing Reynolds stress transport models. Accurate prediction of separated flows are 
important for both external and internal aerodynamics, including off-design conditions, wing-body junction 
flows, turbomachinery diffusers, etc. By isolating key problem areas, we hope to lay the ground work for 
future improvements to RSMs. The long-range goal is the improvement of RANS models for predicting 
separated flows. 


Case no. 

Hump Geometry 

Re(_EisA) 

Simulation 

Reference 

1 

Convergent-divergent channel 

8,404 

DNS 

4 

2 

Periodic hill 

12,714 

LES 

6 

3 

Rounded step 

13,700 

LES 

5 


Table 1. 2-D hump geometries studied using highly-resolved simulations. 


Balance of RSE terms in plane channel flow 



II. 2-D hump datasets 

The hump geometries studied here have flow constrained between two walls, one wall parallel to the 
incoming flow and the other contoured smoothly, to cause flow separation on the leeward side of the hump 
at the given Reynolds number. The 2-D data sets in Table 1 can be found in an online turbulence modeling 
resource. 7 They provide the mean flow, Reynolds stresses, and budgets. Note that dissipation values in the 
LES datasets have been calculated as a remainder of the budget balances (i.e. , to ensure the budget balances 
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are zero), rather than by averaging the instantaneous data, and hence need to be used cautiously. 

The separation bubble sizes for cases 2 and 3 are on the order of the height of the hump. However, case 
1 produces a very small separation with wall-normal height of y + ss 20 and is of shorter length (Fig. 3). 
In the shear layer region just beyond the start of separation, all three cases yield localized V /e far from 
equilibrium. In cases 2 and 3 the maximum localized V /e is near 3-4, whereas in case 1 it is V /e ss 5.5. 



(a) case 1 



(b) case 2 



Figure 3. Geometries of humps studied, with flow streamlines. The U=0 line is shown ( ). 


III. Deficiencies of existing RSMs 

Reynolds stress transport models can be written as: 

T>UjUi „ „ „ 

^ = Pij ~ £ ij T T D Vl ij + D Pt ij + Tij (1) 

As described in the Introduction, the closure models for pressure-strain (<I» tJ ) and dissipation (e^) terms 
in the Reynolds stress equations (RSE) are generally the major sources of error. This is also discussed in 
Bentaleb et al. 2 Although not shown here, the magnitudes of the budget terms in the RSE generally increase 
in the latter half of the separation bubble. Therefore it is instructive to look at Reynolds stress budgets in 
this region. Figure 4 shows the budgets at the X c = 65 station for case 1. There are two regions of large 
magnitude of loss/gain in the stress transport budget; the first in the near-wall region and the second in 
the separated shear layer. It is important that models predict both the near-wall inhomogeneity and the 
non-equilibrium shear layer accurately. 

A. Pressure-strain modeling 

Models for pressure-strain typically assume turbulence is homogeneous. Furthermore, most do not correctly 
account for near- wall effects and hence behave non-physically in near- wall regions. As will be shown below, a 
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Figure 4. Budgets of Reynolds stresses at separation location X c = 65 in the convergent divergent channel 
(case 1). 


priori studies of models using the DNS and LES data sets show models to fail not just near the wall, but also 
near the shear layer region where non-equilibrium is particularly high (V/e > 3). In the study of Bentaleb et 
al., 5 * non-local effects were shown to be pronounced in the separated near-wall region. The elliptic relaxation 
model 8 accounts for the near- wall inhomogeneity of —£ij by an elliptic relaxation of a near- wall correction 
to far-held values. The elliptic correction is solved in a coupled manner and hence cannot be studied a priori; 
full RANS computations have been conducted and are discussed in the last section. 

A priori studies have been conducted using the linear and non-linear pressure-strain models listed in 
Table 2. Homogeneous pressure strain models can be expressed as a non-linear function of stress anisotropy, 
strain and vorticity tensors using the Cayley-Hamilton theorem. A general expression (see Pope 9 ) is: 
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Modelers have considered truncated approximations of this series. The simpler ones in Table 2 are linear 
or quadratic in stress anisotropy and vorticity (Isotropization-of-Production 10 (IP), Launder- Reese-Rodi 10 
(LRR), and Speziale-Sarkar-Gatski 11 (SSG)), and more complex ones are cubic (Shih-Lumley 12 (SL)) and 
quartic (Fu-Launder-Tselepidakis 13 (FLT) and Two-Component-Limit 14 (TCL)). This collection of homoge- 
neous pressure-strain models has been tested on all three of the 2-D separated flows listed in Table 1 through 
a priori analysis. However, we only show results from the rounded step case 3. Results from the other two 
cases support the same conclusions. 
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* Original LRR 10 used C\ = 1.5, but most LRR imple- 
mentations in the literature use 1.8. The < F i j predictions are not significantly influenced by this modification. 


Table 2. Coefficients of 'I', .. Models. 


The values of stress anisotropy, strain, and vorticity from the highly-resolved simulations are used to 
evaluate the pressure-strain models using the algebraic expression of Equation (2). Locally, the boundary 
layer prior to separation accelerates, increasing the strain and vorticity levels. For the rounded step, all the 
models incorrectly predict the 4> 1? components near the wall upstream of separation at X c =-21, as seen in 
Figs. 5a, 6a, and 7a. The cubic model of SL and the quartic model of FLT do not over-respond as much 
as the other models in this region. In the beginning of the separated region, the flow encounters the non- 
equilibrium separated shear layer and increasing inhomogeneity. As seen in Figs. 5b, 6b, and 7b, immediately 
after separation (X c =5), the models are inaccurate away from the wall where P/e is high. However, SL 
predicts the <4>n component in this region more accurately. At stations further downstream well into the 
separation bubble (Figs. 5c and d, 6c and d, and 7c and d), all the models are reasonably accurate in the 
separated shear layer regions ( y/H near 0.8 and 0.6, respectively), but continue to be inaccurate close to the 
wall. The normal components of < f> J j have the wrong sign close to the wall, while the shear component has 
the correct sign. It should be noted that the highest values of pressure-strain are observed at the wall in the 
DNS/LES data sets close to reattachment, rather than in the shear layer. The near- wall maxima/minima 
have been clipped in Figs. 5, 6, and 7 for clarity. 
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The a priori analysis demonstrates that existing pressure-strain models are deficient very near the wall 
as well as in the start of the separated shear-layer region. However, it is difficult to ascertain the particular 
importance of this deficiency, because individual errors from the various modeled budget terms can often 
cancel each other, making it difficult to isolate a primary cause. 15 In any case, as will be shown in Section 
V, a posteriori computations of these 2-D separated flows are poor. 





Figure 5. Pressure-strain component <3>n predicted by a priori analysis of various models on the rounded step 
(case 3) geometry. The maximum values of V/e at each of these streamwise locations are given. 


B. Dissipation modeling 

From the budgets in the separated region of the convergent-divergent channel (Fig. 4), it is clear that 
dissipation plays an important role relative to other budget terms near the wall. The normal components 
£n and £33 are large and are an important part of the balance in near-wall turbulence. Plotting the trace of 
dissipation ( Eh/2 ) in Fig. 8, the dissipation is seen to rapidly increase half-way into the separation bubble, 
mostly near the wall. This growth of dissipation likely influences flow behavior near reattachment. 

Many near-wall dissipation models have been developed. Here the anisotropy models of Lai and So 16 
(LS) and Hanjalic et al. 1 ' (HJH) were tested. These models use different blending functions sensitized by 
the wall normal vector to blend the isotropic far field to the anisotropic near-wall dissipation. The models 
are described in Equations (3) and (4). 

LS 16 : 
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Figure 6. Pressure-strain component <£>12 predicted by a priori analysis of various models on the rounded step 
(case 3) geometry. 






Figure 7. Pressure-strain component $22 predicted by a priori analysis of various models on rounded step 
(case 3) geometry. 
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The HJH model predicts Eij somewhat better than LS both near and away from the wall as seen in the 
plane channel flow computation in Fig. 9. The isotropic dissipation curve is also included for comparison, as 
isotropy is a widely used assumption. 



Figure 8. Dissipation profiles (ea/2) at various separa- 
tion coordinates X c locations from DNS of convergent- 
divergent channel (case 1). A rapid increase in near- 
wall dissipation is seen in the latter half of the separa- 
tion bubble. 


Figure 9. Dissipation tensor predicted using the 

anisotropic models of HJH (- - -) and LS ( • ) 

in a plane channel flow at i?e r =590 are compared with 
DNS ( ). 


IV. Explicit algebraic Reynolds stress models (EARSM) 

EARSMs are computationally inexpensive and typically more robust compared to Reynolds stress trans- 
port models. They retain the modeling components of the RSE, without the dynamics of stress transport. 
The modeled transport equation for the anisotropy tensor can be written as: 

UiUj UjUj „ 

ij k~ D ) = k~ ^ ~ + Pij ~ £ij + ® ij 

where is the combined effect of turbulent transport and viscous diffusion. 18 An implicit algebraic stress 
relation is obtained by assuming that the terms on the left hand side are zero, and an explicit algebraic 
Reynolds stress relation can be found by expressing the anisotropy tensor in terms of specified tensor bases. 

In the Gatski and Rumsey version of EARSM, 18 a linearized SSG 11 pressure-strain model is used with 
a different assumption for the diffusion term D ij. In the Wallin and Johansson version of EARSM, 19 a 
modified version of the LRR 10 pressure-strain model is used. Both models assume dissipation to be isotropic 
= |ei 5ij. Near the wall the viscous and turbulence dissipation are the dominant balances, leading to 
the limiting behavior of Eij — > UiUj£/k , i.e. the anisotropies of dissipation and Reynolds stresses are the 
same. However, far from the wall the dissipation tensor is isotropic. The assumption of isotropy is fairly 
accurate in most regions of flows, except near the high V / e region of separated shear layers. The Wallin and 
Johansson 19 model was also modified and tested using the assumption of dij = ciij. 

A priori testing indicates that the three models under-predict anisotropy in the near-wall region prior to 
flow separation, as seen in Figs. 10 and 11. As noted earlier, the homogeneous models are inadequate 
to predict the near-wall damping of the normal stress. After separation the two EARSMs 18, 19 predict the 
trends in anisotropy (away from the wall) reasonably well. However, the anisotropic dissipation EARSM 
(with dij = ciij) does not bring an improvement in near- wall anisotropy compared to Wallin and Johansson’s 
original EARSM. 
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Figure 10. Anisotropy of shear-stress at different streamwise coordinates in the separated region of the rounded 

step (LES ) predicted by EARSM of Gatski and Rumsey 18 ( — ), Wallin and Johansson 19 ( ) and 

its variant that assumes dissipation anisotropy to be same as that of the Reynolds stress (■■ ■ ■). 






Figure 11. Anisotropy of principal normal-stress at different streamwise coordinates in the separated region 
of the rounded step. Line types are those of Fig. 10. 
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V. Full computations 


The computations for all the cases were conducted using OpenFOAM®. 20 The mean and turbulent 
equations are solved using the bi-conjugate gradient linear solver, with residuals converged to eight orders 
of magnitude. A blended upwind and central differencing was used for discretization of advection terms. 
To investigate the effect of accounting for non-local effects close to the wall, the elliptic blending model of 
Durbin 8,21 was tested. It had to be solved using an unsteady formulation for the solution to be tractable. 
Even so, it only converged adequately for case 3, but not for case 1. In this implementation, the elliptic 
equation is solved decoupled from the stress equations. 

All the models tested here use the isotropic dissipation tensor and the LRR pressure-strain model. As 
all pressure-strain models in the a priori tests of Section IIIA yielded similarly poor results near the start 
of the separated shear layer, LRR can be considered representative for the purposes of this testing. Two 
models were used to account for near-wall effects of inhomogeneity: elliptic relaxation 8,21 and Launder- 
Gibson, 22 the latter of which is sensitized using wall-normals. Both of these RSMs solve the dissipation (e) 
transport equation along with an e wall-function at the boundaries. Pressure-diffusion, D p ij, and turbulent 
diffusion, X) ? - , are treated along with the viscous diffusion using the generalized gradient hypothesis of Daly 
and Harlow. 23 The quality of the results are assessed primarily by the ability to predict the separation and 
reattachment locations. 

First, for case 1, the convergent-divergent channel, the grid used for the 2-D RANS computations (288 x 
193) is 8 times coarser along the streamwise index and twice coarser along the wall normal direction than 
the original DNS grid (2304 x 385). The wall-normal mesh has a y + ss 1 at the wall based on inlet u T . The 
DNS predicts a separation bubble less than one unit in length, with separation at Xg ~ 5.8. As seen in 
Fig. 12 and Table 3, both EARSM 19 and Launder-Gibson RSM 22 predict delayed separation and too large 
a bubble extent. 

The 2-D rounded step case 3 is solved on the same grid used by the LES 5 (769 x 161). This grid was 
designed to yield wall spacing of y + < 1. The models all over-predict separation bubble extent, as seen in 
Fig. 13 and Table 3. EARSM predicts similar separation location and topology, but the bubble extent is 
significantly larger than the reference LES. The Launder-Gibson RSM shows a reattachment that is even 
further delayed with the separation start Xg also predicted significantly too far downstream. An abnormal 
reattachment behavior is also observed, where the streamwise flow reverses twice along the wall-normal 
direction. This behavior is known to be caused by the modeling of the e transport equation. 2 The elliptic 
relaxation RSM predicts the separation start somewhat late and also overpredicts separation extent, in spite 
of accounting for near-wall inhomogeneity in <Fy — by solving an additional elliptic equation for each 
stress component. 

Example wall skin friction coefficient plots for the two cases are shown in Fig. 14. Here, only EARSM 19 
results are compared to the reference solutions in order to give an indication of typical RANS behavior for 
these cases. For case 1, the RANS predicts separation too late, along with a separation bubble that is too 
long. For case 3, the RANS predicts separation location well but the separation bubble is again too long 
and reattaches too late. Both the convergent-divergent channel and the rounded step cases have proven 
to be challenging for today’s models. Although not shown here, widely-used two equation models such as 
Menter’s Shear Stress Transport (SST) 24 also over-predict separation bubble extent. The periodic hill is 
similarly problematic for RANS models. In terms of accurate separation bubble prediction, isolating and 
treating only one of the budget models accurately — such as pressure-strain or dissipation tensor — is likely 
to be insufficient. 


Laval Hump (case 1) Rounded Step (case 3) 




X R 


Xr 

DNS/LES 

5.76 

6.60 

0.83 

4.35 

EARSM 

6.31 

8.61 

0.91 

5.00 

Launder Gibson RSM 

7.21 

8.37 

1.50 

5.35 

Durbin RSM 

- 

- 

1.09 

4.82 


Table 3. Separation and reattachment locations for case 1 and 3, non-dimensionalized by half-channel and 
hump height respectively. 
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Figure 12. Separation streamlines and flow reversal line ( ) predicted by DNS and RANS models for 

the convergent-divergent channel (case 1). 



Figure 13. Separation streamlines and flow reversal line ( ■ ) predicted by LES and RANS models for 

the rounded step (case 3). 
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(a) Convergent-divergent channel (case 1) (b) Rounded step (case 3) 

Figure 14. Skin friction coefficient predicted by the EARSM 19 for two hump cases. 


VI. Conclusion 

Using fully-resolved simulation data for 2-D separated flows, a priori studies of existing pressure-strain 
models show that all models fail both near the wall and near the start of the separation bubble. There is 
a possible correlation between the poor predictive capability of the pressure-strain models in the high non- 
equilibrium shear region near the start of the bubble and the poor predictions of the full RANS computations 
for separated flows. However, there are many other modeling factors that are potentially influential as well, 
which are difficult to isolate. Further testing of dissipation models is necessary to ascertain their strengths 
and weaknesses for 2-D separated flows. Dissipation increases rapidly towards the reattachment region of 
the flow; hence dissipation models likely play an important role in predicting the flow physics of the bubble. 
Accurate near-wall modeling will also likely prove to be critical. This research represents a first step toward 
identifying problems with existing pressure-strain and dissipation models applied to separated flows. 
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